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ABSTRACT 

We study the short-term effects of an interniediate mass black hole (IBH) on the 
orbit of star S2 (S02), the shortest-period star known to orbit the supermassive black 
hole (MBH) in the centre of the Milky Way. Near-infrared imaging and spectroscopic 
observations allow an accurate determination of the orbit of the star. Given S2's short 
orbital period and large eccentricity, general relativity (GR) needs to be taken into 
account, and its effects are potentially measurable with current technology. We show 
that perturbations due to an IBH in orbit around the MBH can produce a shift in 
the apoapsis of S2 that is as large or even larger than the GR shift. An IBH will also 
induce changes in the plane of S2's orbit at a level as large as one degree per period. 
We apply observational orbital fitting techniques to simulations of the S-cluster in the 
presence of an IBH and find that an IBH more massive than about 1000 M0 at the 
distance of the S-stars will be detectable at the next periapse passage of S2, which 
will occur in 2018. 

Key words: black hole physics ~ stellar dynamics - methods: TV-body simulations - 
Galaxy: centre 



1 INTRODUCTION 

The innermost arcsecond of the Milky Way harbours a clus- 
ter of young massive stars (the S-star cluster) in eccentric 
orbits around the supermassive black hole (MBH). Near- 
infrared observations of the cluster allow a precise deter- 
mination of the trajectories of about 20 stars. These can 
be used to derive fundamental parameters like the mass 
of the MBH and the distance to the Galactic centre, as 
well as to constrain the gravitational potential and test 
predi c tions from general r elativity (GR) (|Rubilar fc Eckartl 
I2OOII : IZucker et all I2OO6I ). Since the collapse of a giant 
molecular cloud at the distance of the S-cluster is pre- 
vented by the tidal force of the MBH, it is believed 
that either the S-stars formed in-situ via a non-standard 
process or that they originated outside the black hole's 
sphere of influence and then migrated to their current lo- 
cation. Among the in-s itu models are the two disks model 



Lockmann et al.l 200!^) and the eccentric instability model 



Madigan et al 



I2OO9). The two main models suggested 
for the transport of the stars are the cluster infall sce- 
nario ifCerha rd 200 li). also aided by an interm e diate-mass 



black hole flBH) (|Hansen fc MilosavlievidlioO^ : iKim et ai] 



2004 ; iLevin et al.]|2005l : iMerritt et al.ll2009l: jFuiii et al.ll2009l. 



201c ) , and the binary disrupt ion scenario (|Gould fc QuiUenI 



20031: iPerets et al.ll2007l . l2009l ). The properties and efficienc y 
of these models are discussed in lPerets fc GualandrisI l|2010l ). 

The star with the shortest orbital perio d (~ 15 yr), 
called S2 (jSchodel et al.ll2002l : [Chez et al.ll2003h . has a semi- 
major axis a = (0 .1246 ± 0.0019)" an d an eccentricity 
e = 0.8831 ± 0.0034 l|GiUessen et al.ll2009l ). Adopting a dis- 
tance to the Galactic centre _Ro = 8.28 kpc, a ~ 5.0018 mpc, 
while the periapsis and apoapsis distances are, respectively, 
rp = 0.585 mpc, Va = 9.419 mpc . For this set of orbital 
elements, GR precession is measurable with current instru- 
ment ation on a time-scale of about 10 years (|Gillessen et al.l 
I2OO9I ). 

The relativistic (prograde) advance of the periapse an- 
gle is given in the case of a non-rotating black hole by 
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per radial period (|Weinberel Il972t ). where -Rs — 
2GMmbh/c'^ is the Schwarzschild radius of the black hole. 
For the Galactic centre black hole, Rs ~ 4. 1 x 10~^ pc for 
an as sumed mass of Mmbh = 4.3 x 10^ Mq (|Gillessen et al.l 
[iooi). As a result of the orbit's precession, there is a dis- 
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placement in the star's apoapse position that is given by 

67rGA/MBH 



Ara ^ a(l + e)Azu : 



c2(l - e) 



(2) 



again per radial period; note that this expression is indepen- 
dent of the semi-major axis. As seen from Earth, this shift 
corresponds to an angle on the sky of 



0.097 mas 



1 



4.3 X 106 Mq 



..28 kpc 



(3) 

which amounts to ~ 0.83 mas for star S2. 

A complicating factor is the likely presence of a dis- 
tributed mass within S2's orbit, consisting of stars and stel- 
lar remnants. Assuming that the mass density follows r 
with r the distance from the MBH, the advance, per radial 
period, of orbital periapse due to the distributed mass is 



2tv 



jVA(a) 
Mmbh 
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(4) 



in the retrograde sense; here M*(r) is the distributed mass 
enclosed within rad ius r and F = (3/2, 1) for 7 = (0, 1) (e.g. 
iMerritt et allbOlOl '). Setting F x 1, this implies, for S2, a 
shift on the sky of 



AQa ~ 0.69mas 



Mmbh 



(5) 



This is comparable with the relativistic precession if the dis- 
tributed mass within S2's orbit is ^ 10~"^AfMBH. Disentan- 
gling these two sources of precession will be difficult, requir- 
ing measure ments of the radial ve locity of S2 near periapse 
passage (e.g. lAneelil fc Sahall2010h . 

Non-spherically-symmetric perturbations, if present, 
can affect not just the periapse angle but also the angu- 
lar momentum of S2's orbit, resulting in changes in e and 
in the orientation of its orbital plane. Potentially the largest 
source of such perturbations is a second massive black hole 
orbiting somewhere near S2. In fact, the orbits of the S- 
stars are consistent with the long-term presence of an IBH 
in their midst, with mass ~ 10^ (jMerritt et al.l l2009l : 
iGualandris fc MerrittllioogI ). While an IBH would also con- 
tribute to the periapse advance of S2, its key signature would 
be a change in S2's orbital angular momentum. 

In this work, we combine high-accuracy A''-body sim- 
ulations (§2) with orbital fitting techniques (§3) to investi- 
gate the observable effects of an IBH on the orbit of star 
S2 over a time-scale of a few orbital revolutions. In §4 we 
discuss the other types of perturbation that could poten- 
tially induce changes in S2's orbital angular momentum and 
compare them with perturbations from an IBH. §5 sums up. 



which appears to be unbound. We are left with a sample of 
21 stars with well defined orbits, for which we determined 
positions and velocities at 2008 AD from the classical ele- 
ments. The masses of the S-stars were set to 10 M0 (e.g. 
lEisenhauer et aLlbOOSl ) except for star S2 f or which a value 
of 20 Mq was adopted ()Martins et al.ll2008h . 

The IBH is placed on a Keplerian orbit around the 
MBH. We adopt four different values for the mass ratio of 
the black hole binary q = Mmbh/Mibh = (1-0 x 10"*, 2.5 x 
10"^, 5.0 X lO""*, 1.0 X 10"^), five values for the semi-major 
axis a = (0.3, 1, 3, 10, 30) mpc, four values for the eccentric- 
ity e = (0,0.5,0.7,0.9) and twelve choices for the direction 
of the orbit's angular momen tum vector (the same set as in 
IGualandris fc MerrittI ()2009l l). for a total of 960 sets of ini- 
tial conditions. The IBH begins from orbital periapsis in all 
cases. The initial value of the mean anomaly is likely to be 
unimportant in all cases for which the orbital period of the 
binary is much smaller than the integration time (50 yr), i.e. 
a = 0.3, 1, 3 mpc. For a = 30 mpc, only part of the IBH orbit 
is sampled by the integration, and in principle the choice of 
the initial position might have an effect on the interaction 
with S2. However, at a distance of 30 mpc the IBH is com- 
pletely outside S2's orbit for e ;£ 0.7, and there is essentially 
no detectable signature, as shown below. The only set of 
simulations for which the initial mean anomaly of the IBH 
might be important is the a = 10 mpc, which corresponds to 
an orbital period of about 43 yr for the IBH. We perform an 
extra set of simulations for a = 10 mpc starting the IBH at 
apoapsis rather than periapsis and we compare the results 
in the two cases. 

We advanced each A^-body system in time using the 
AR-CHAIN code l|Mikkola fc Merrittlliooi ) . a recent imple- 
mentation of the algorithmic regularisation method that is 
able to reproduce the motion of tight binaries for long pe- 
riods of time with extremely high precision. The code com- 
bines the use of the chain structure, introduced originally by 
iMikkola fc AarsethI (|l993 l). with a new time transformation 
to avoid singularities and achieve high precision for arbi- 
trary mass ratios. Note that we self-consistently follow not 
just the interactions of S2 with the MBH and IBH, but all 
other interactions as well, including star-star interactions. 
The integration interval was 50 years, during which time S2 
performs three full orbits. 

The AR-CHAIN code includes relativistic corrections 
to the accelerations up to 2.5 post-Newtonian order for all 
interactions involving the MBH particle. General relativistic 
advance of the periapse, which operates on a time-scale 



2 INITIAL MODELS AND NUMERICAL 
METHODS 

We consider A^-body models that include a MBH, an IBH 
and the S-star cluster. The initial conditions for the stars are 
derive d from the orbital elements given by I Gillessen et al.l 
||2009D. For star S2, we take the improved elements from 
iGiUessen et al] l|2009l ). Of the 28 stars for which those au- 
thors provide orbital elements, we exclude the six stars (S66, 
S67, S83, S87, S96 S97| which hkely belong t o the clockwise 
disk l|Genzel et al.ll2003l : |Paumard et aLllioO^ and star Sill 
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is accounted for by the IPN and 2PN terms. The dissipa- 
tive term arising from the emission of gravitational waves 
is accounted for by the 2.5PN term; this term is potentially 
important for the IBH, for which the associated coalescence 
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time-scale is 
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is the reduced mass of the IBH/MBH binary. The orbital de- 
cay time-scale for the black hole binary is shown in Figure [1] 
for the two extreme values of the mass ratio and the four 
adopted values of the initial eccentricity. This time-scale is 
always much longer than our integration interval of ~ 50 
yr. In addition, it is longer than the main-sequence lifetime 
(~ 10^ yr) of a 20 solar mass star for all initial configura- 
tions excepting the cases a — 0.3 mpc and e ^ 0.7. In the 
former runs, it is justified to associate our initial parame- 
ters for the IBH/MBH binary with the parameters at some 
much earlier time, e.g. the epoch preceding formation of the 
S-stars. In the latter runs, the orbit of the IBH at some much 
earlier time would have been larger and/or more eccentric. 
The maximum relative variation of the binary semi-major 



axis in the A'^-body integrations is Aa/a ~ 
absolute variation of the eccentricity is Ae > 



3 IBH PERTURBATIONS 
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In the Schwarzschild metric, the argument of periapse zu 
evolves due to in-plane precession. The two remaining angles 
that define the orientation of the orbit, i, the inclination, 
and n, the position angle of the ascending node, are fully 
conserved in the relativistic two-body problemQ The semi- 
major axis and eccentricity are conserved at the IPN level, 
and we expect very small devia tions due to higher order FN 
corrections in the integrations (|Soffellll989l ). In the limit of 
small star-to-black hole mass ratio, the semi-major ax is and 
eccentricity in the FN approximation are given by (ISoffell 
1 19891 ) 
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We follow the standard practise of using the plane of the sky 
as the reference plane for defining (f2,i). 
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Figure 1. Time-scales associated with orbital evolution in our 
models. Solid lines show the GW time-scale, Eq. ((Tjl, for a black 
hole binary with q = 10"* (top) and q = 10~^ (bottom), for four 
different values of the eccentricity e = 0,0.5,0.7,0.9. Dashed- 
dotted lines show the GR precession time-scale for two difi'erent 
values of the eccentricity: e = (upper line) and e = 0.9 (lower 
line). The vertical dotted lines represent the adopted values for 
the binary initial semi-major axis. The filled grey region indicates 
the estimated ages of the S-stars while the striped area shows the 
radial range of S2's orbit. 



is the specific post-Newtonian energy and 

3GM 



1 + 
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the specific angular momentum. Here, and are the 
relative position and velocity vectors between the star and 
the MBH, M is the total mass, and c is the speed of light. 

The presence of the other S-stars introduces small devi- 
ations from spherical symmetry in the gravitational poten- 
tial but the effect on the orbital elements over these short 
time-scales is negligible, as we show in Section [5] There- 
fore, the variations that we observe in the orbital elements 
of star S2 in our simulations, namely the semi-major axis, 
eccentricity, inclination, and position angle of the ascending 
node, can be attributed to perturbations from the IBH. 

Figure[2]summarises the changes in the orbital elements 
of S2 found in the A'^-body integrations. Flotted are the vari- 
ations over one revolution averaged over the twelve different 
orientations of the initial IBH/MBH orbit. The dotted lines 
in the periapsis and apoapsis panels indicate the variations 
due to GR. The shift in the periapse corresponds to an ob- 
servable angle AOp — AOa (1 — e) / (1 -f- e), where AOa is 
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Figure 2. Average changes in tlie orbital elements (semi-major 
axis, eccentricity, inclination, position angle of the ascending 
node, periapsis and apoapsis) of star S2 over one full orbit, versus 
the mass ratio of the black hole binary. Different symbols are for 
different initial semi-major axes of the binary. Each point is an 
average over the 12 orientations of the IBH/MBH orbital angular 
momentum vector. The dotted lines represent the GR shift in the 
periapse and apoapse. 



defined in Eq. [S] Note that the variations in the inclination 
i and position angle il of the ascending node reach values 
close to 1 degree for the most massive IBHs considered. This 
is of the same order as the current observational accuracy 
(~ 0.7 deg) . In 50 years, this value will drop to about 0.4 deg, 
assuming there are no technological improvements. 

While precession induced by the PN terms is restricted 
to the orbital plane, an IBH induces more general changes 
in the orbital elements, including changes in the direction 
of the orbital angular momentum vector. We measure the 
latter via the angle 
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Figure [3] plots for all the runs, after averaging over the 
12 different IBH orientations. 

The changes in the orbital plane of S2 are larger for 
more massive IBHs and can reach values of ~ 1 degree for 
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Figure 3. Mean variation of the orbital plane of star S2 as a 
function of the binary mass ratio. For each combination of bi- 
nary mass ratio, semi-major axis and eccentricity, the results are 
averaged over the 12 orbital orientations of the IBH. 



of the MBH/IBH binary orbit such that changes are largest 
for a ^ 10 mpc. 

The three-body system compo sed of MBH , IBH and 
S2 is reminiscent of a Kozai triple l|Kozailll962l ). However, 
the changes that we observe in the orbital elements are not 
generally attributable to the Kozai mechanism. Kozai oscil- 
lations can be induced if (i) the MBH - S2 - IBH system 
can be regarded as a hierarchical triple, with well-separated 
orbital periods for the inner and outer orbit; (ii) the period 
predicted for the Kozai oscillations is shorter than that of 
any other precessional period, in this case, GR precession; 
(iii) the outer orbit is largely inclined with respect to the 
inner orbit. The first condition is only satisfied in the runs 
with a — 10, 30 mpc. 

The timescale for Ko zai oscillations can be written as 
(jKinoshita fc Nakai|[2007l ) 
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Out-of-plane motion is also affected by the size 



3\/67r Finn MlBH 

where Pinn and Pout are the period of the inner and outer 
binary, eout is the eccentricity of the outer orbit and IC is 
a numerical coefficient which depends only on the initial 
values of the relative inclination angle a, the inner binary 
eccentricity einn and the inner binary argument of periapsis 
■ou. The maximum eccentricity emax attained by the inner 
binary is also a function of the initial values of a, einn and 
zu only. Figure [3] shows Tk as a function of the black hole 
binary mass ratio for the assumed values of the eccentricity 
and for two allowed separations. We find that Tk is always 
longer then S2's GR precession timescale of 2.5 x 10^ yr for 
a = 30 mpc. For the remaining case a = 10 mpc, Tk is short 
enough only for q ^ 2.5 x lO"** and e ^ 0.7. This restricts the 
applicability of the Kozai mechanism to a small subset of our 
simulations, in contrast with the results shown in Figure [2] 
and in the following section. We conclude that the Kozai 
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Figure 4. Timescale for Kozai oscillations in star S2 as a function 
of the black hole binary mass ratio for different values of the 
binary eccentricity. The left panel refers to simulations with a = 
10 mpc while the right panel is for a = 30 mpc. The horizontal 
lines marks the GR precession timescale for S2. 



mechanism is not the dominant effect producing changes in 
orbital e and i in our simulations. In those cases where it is 
potentially relevant, the eccentricity is predicted to oscillate 
between Cmin = 0.24 - 0.74 and Cmax = 0.89 - 0.99 over 
one Kozai period, depending on the values of a and vj. This 
implies variations of the order of 2 x 10~^ — 6 x 10""^ over 
the integration time of 50 years (for a — 10 mpc). 

The apoapsis shift due to perturbations from the IBH 
is very sensitive to the binary parameters. In the case of 
g ^ 5 X 10~* and a £ 3 mpc, the shift over one revolution due 
to the IBH becomes larger than the relativistic shift. This 
suggests a variation in the orbital elements which is poten- 
tially observable with current instrumentation. 

However, the observability of variations in the orbit of 
S2 depends on several factors. In the following section we 
thoroughly examine all such factors and use orbital fitting to 
determine whether an IBH is detectable via on-going mon- 
itoring of the S-cluster. Theoretically, it would be possible 
to use other S-stars to investigate the effects of a hypothet- 
ical IBH. Given that the shifts in the apparent location of 
periapsis and apoapsis depend only on the eccentricity and 
not on the semi-major axis of the stellar orbit, it would 
seem appropriate to consider all stars with e^O.8 for such 
an analysis. From an observational point of view, however, 
S2 is the only star in the sample which is bright enough and 
not affected by confusion to allow for meaningful tests of the 
gravitational potential. We therefore limit our study to star 
S2. 



4 ORBITAL FITTING 

In this section we extract observational-like data from the 
simulated orbital traces of S2. We assume that eight or nine 
astrometric epochs can be obtained each year over the course 
of 50 years. The eight or nine yearly epochs are not evenly 
distributed but are spread over seven months only, thus tak- 
ing into account the fact that the Galactic Centre is acces- 
sible with NIR observations only for part of the year. At 
the chosen epochs, the original, simulated positions are per- 
turbed by an astrometric error, assumed to be distributed 
in a Gaussian fashion. We use a value of 300 fias per coordi- 
nate, which is a conservative assumption for S2 (|Fritz et al.l 
I2OI0I ) . The statistical uncertainty of the measured S2 posi- 



tions is smaller than the assumed value. However, unrecog- 
nised confusion events with fainter stars are an additional, 
important error source, such that we consider our value re- 
alistic. Furthermore, we assume that the radial velocity of 
S2 is determined at two epochs per year. Here, we adopt a 
Gaussian error of 15kms~^. This value is typically reached 
with current NIR medium-resolution spectrographs. 

In this way we obtain 960 simulated data sets that 
are reasonably close to what one would obtain by simply 
continuing the monitoring of orbits in the Galactic Centre 
with existing instruments. The fact that within the next 50 
years new facilities with improved angular resolution will 
become available (such as the NIR interferometers ASTRA 
(|Pott et al. II2OO8I ) and GRAVITY l|Eisenhauer et al.ll2005l ). 
or the extremely large telescopes TMT and E-ELT) means 
that an IBH probably will be detectable more easily than 
what we derive here. There is a second reason why we con- 
sider our procedure conservative. Observers may adopt a 
sampling strategy for the orbits other than a simple con- 
stant rate, as we assume here. Given the fact that the most 
constraining part of an orbit is the periapse passage, and the 
event is predictable, an intensification of the observations 
around the periapse passage provides an improved sampling 
pattern. 

We fit each of the 960 da ta sets with the same code as 
used in lCillessen et aL 

I (|2009l). From the fits, we determine 
the full set of 13 parameters describing the orbits and the 
potential: The six orbital elements (a, e, i, Q, vj, tp) and 
seven parameters describing the gravitational potential of 
the MBH: Mass, distance, on-sky position (2 parameters) 
and velocity (3 parameters) of the MBH. While in the simu- 
lations these quantities are known, this is not the case when 
the data sets are considered as mock observations. Then 
these quantities have to be treated as free fit parameters, 
since they need to be determined from the same orbital data 
from which the presence of an IBH shall be judged. Cur- 
rently, nearly all constraints on these parameters actually 
come from the S2-orbit which we consider here. In future, 
some of the parameters describing the MBH might be deter- 
mined independently from the orbit of other S-stars. We ne- 
glect this here and keep all seven parameters completely free. 
This is conservative, because additional constraints would 
make the presence of an IBH more easily detectable. 

Using a purely Keplerian point-mass model is inade- 
quate for all our data sets. The relativistic precession is 
too large during the 50 years of evolution and therefore we 
include the first order PN corr ection to the e quations of 
motion when fitting the orbits (|Gillessen et al.l 2009). For 
example, all twelve data sets with {q = 1.0 x 10~^, a = 
30 mpc, e = 0.9) are well fit by the relativistic equations. 
In contrast, none of these data sets can be described by a 
purely Keplerian point-mass potential. Still, the relativistic 
potential only yields a perfect fit for some of the data sets. 
For others, the IBH perturbs the dynamics too strongly, and 
would therefore be detectable. Given that our sampling in 
radial velocity is rather sparse, the inclusion of special rela- 
tivistic corrections to the radial velocities is not needed. 

Note that it is not legitimate to assess the goodness of a 
fit by comparing the obtained parameters with those used as 
input for the simulations. Instead, we use the reduced for 
each of the 960 fits to decide whether a fit is acceptable or 
not. We obtain values between 0.88 and 357. In addition, we 
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Figure 5. Examples of fit residuals. Top row: A fit using a sim- 
ulated data set with (q = 1.0 X 10"*, a = 3mpc, e = 0). Middle 
row: a fit witli (q = 2.5 X 10"*, a = 0.3 mpc, e = 0.9). Bottom 
row: a fit with (g = 5 X 10"*, a = 1 mpc, e = 0). The first (top) 
one is classified as acceptable and has a reduced = 1.15. The 
other two are classified as not acceptable and have reduced 
values of 4.4 and 22.7 respectively. 



examine the residuals of each fit by eye, dividing them into 
two categories: A set of fits for which the residuals do not 
visually show obvious correlations and a set for which it is 
apparent that the chosen gravitational potential model is not 
adequate. Figure [5] shows typical examples of the residuals 
from three simulations. 

Each of the simulated data sets has 418 ± 1 astro- 
metric data points and 142 ± 1 radial velocity measure- 
ments. The exact numbers vary due to the random vari- 
ations of the sampling pattern. Given that our fits have 
13 free parameters, the number of degrees of freedom is 
ridof « 2 X 418 + 142 - 13 = 965. 

Figure [E] shows the distribution of reduced together 
with the flag whether a given fit is acceptable or not. Clearly, 
the reduced can be used as discriminator. The largest re- 
duced corresponding to a fit classified as 'acceptable' is 
1.36, the smallest reduced corresponding to a fit classified 
as 'not acceptable' is 1.12. The optimum cut is at a reduced 
of 1.22, yielding as many good fits above as bad fits be- 
low the threshold. The total number of fits misclassified by 
this cut is 41. The total number of bad fits is 409 and corre- 
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Figure 6. Distribution of reduced x ^or the 960 fits in loga- 
rithmic bins. Entries are coloured according to our visual clas- 
sification of whether the residuals of any given fit is acceptable 
(blue/solid) or not (red/hatched). The black, dashed line marks 
the optimum cut at 1.22 separating good from bad fits by the 
value of their reduced x^ . 



spondingly 551 fits have a reduced x^ below the threshold. 
Hence, the IBH would be detectable from the data in ~43% 
of the cases. We find a dependence of the fraction of bad fits 
on the assumed mass for the IBH. For q = lO""*, 2.5 x 10~^ 
5 X 10~^, 10~^ the percentage of detectable IBHs is 15%, 
39%, 51%, and 66%, respectively. 

Having a reduced x^ ^ 1-22 corresponds formally to 
a 4.7cr detection of the effects of the IBH for ridof ~ 965. 
In reality, there will be perturbing effects such as confu- 
sion events, recognised or unrecognised. Hence, the actual 
X^ might be larger than in the simulations and the simple 
Gaussian statistics will not apply. We are nevertheless con- 
fident that by observing the astrometric (and photometric) 
residuals as a function of time the effects of an IBH can 
be disentangled from confusion events. The patterns of the 
residuals contain more information than a single x^ value. 

In Figure [7] we examine the dependence of the reduced 
X^ on the black hole binary parameters. The left panel in 
the first row shows that the fits on average get worse when 
the mass of the assumed IBH is increased, as expected. The 
right panel in the first row shows that a value for the semi 
major axis of the IBH around ~ 1 mpc leads on average to 
the worst fits. These two parameters correlate actually best 
with the reduced x^, and in Figure|S]we show the reduced x^ 
in the MiBH-a-plane. The unacceptable fits occupy a well- 
defined region in this plot. 

The goodness of the fits, on the other hand, is fairly in- 
dependent of the eccentricity of the IBH orbit (second row 
of Figure [T] left panel) . The right panel in the second row 
shows that the periapsis distance p = a (1 — e) of the IBH is a 
less good predictor for the fit than the semi major axis. The 
reduced x^ does not correlate with the finite set of orbital 
orientations probed, as shown in the third row, left panel of 
Figure [7] Here, a represents the angle between the angular 
momentum vectors of S2 and the IBH at the start of the sim- 
ulations. This incidentally argues against the possibility that 
the observed deviations are due to Kozai oscillations, since 
the mechanism requires large relative inclinations to oper- 
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Figure 7. Distributions of reduced for tlie 960 fits as a func- 
tion of the black liole binary initial parameters: mass, semi-major 
axis, eccentricity, periapse distance, orbital orientation and min- 
imum distance to S2. 
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Figure 8. Reduced for the 960 fits as a function of mass 
and semi major axis of the IBH. The plot shows the median at 
each grid point of the values. The thick dashed line marks our 
threshold of x^ = 1-22. Fits right of the line are not acceptable 
and thus the IBH would be discoverable. 
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Figure 9. Reduced x^ as a function of black hole mass (left) and 
eccentricity (right) for two set of runs with the IBH starting at 
periapsis (dots) and apoapsis (circles). 



ate. However, fits with the same sets of parameters {q, e, a) 
but different a can have different values of x^- This means 
that the knowledge of {q, e, a) is insufficient to predict the 
reduced but information on the sky position is necessary. 
Finally, the right panel in the third row of Figure [7] shows 
that the minimum 3D distance between S2 and the IBH also 
correlates with the reduced x^- In general it holds that the 
smaller the minimum distance, the worse the corresponding 
fit. Clearly, this parameter is not independent of the semi- 
major axis. 

The initial parameters adopted for the black hole binary 
are not sampled homogeneously. This is obvious for the first 
three panels in Figure[7]showing the reduced x^ as a function 
of the binary parameters {q, e, a). But it also holds (and is 
less obvious) for the plot investigating the minimum distance 
between S2 and the IBH. 

A comparison of the goodness of fit for the runs starting 
with the IBH at periapsis and apoapsis is shown in Figure [Q] 
In both cases, the semi-major axis of the black hole orbit is 
a — 10 mpc. We find a modest worsening of the in the 
case of an IBH initially at the apoapsis of its orbit. This can 
be attributed to the fact that S2's apoapsis, where the star 
spends most of its time, is about 10 mpc. 

Finally, we also investigated the minimum time required 
for the IBH to become detectable. For this purpose, we re- 
peated the orbital fits for a few cases assuming that the 
observations span 10,15, 20,25, 30, 35, 40, 45 or 50 years 
(Figure llOfl . Our initial conditions are such that the first 
periapse passage of S2 happens after 10 years, the second 
after 26 years and the third after 42 years. Figure [TO] shows 
that the reduced x^ starts to increase beyond our threshold 
of 1.22 after the second periapsis passage for fits that show 
a large reduced x^ after 50 simulated years (red/dashed, 
blue/long-dashed and green/solid curves). Only for fits that 
after 50 simulated years have a reduced ~ 3 is the thresh- 
old passed after the third periapse passage (black curve). 
The discrete nature of periapse passages also is the reason 
why the blue/long-dashed and green/solid curves level out 
after about 45 years. 

This is interesting in comparison with the current status 
of the observations. Since by now only one periapse passage 
of S2 has been observed in 2002, one would not expect to 
have detected an IBH from the actual S2 data so far. This 
is particularly true since the assumed level of accuracy was 
not reached during the first years of the observations (1992 - 
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Figure 10. Reduced a function of tlie lengtli of tlie obser- 

vations for four cases: blue (long-dashed): g = 5 X 10~*, e = 0.9, 
a = 1 mpc; red (dashed): q = 10~^, e = 0.9, a = 1 mpc; green 
(solid): q = 10~^, e = 0, a = 0.3 mpc; black (dot-dashed): 
q = 10^'^, e = 0.7, a = 3mpc. The vertical lines indicate the 
second and third periapse passage covered by the simulations, 
the first one happens at t ^ 10 yr. The horizontal, dashed line is 
our optimum cut at 1.22. 



2002) , and radial velocity information is only available after 
2002 (with the exception of one point in 2000). Hence, the 
first real chance to detect an IBH will be after the next S2 
periapse passage, which will happen in 2018. 



5 OTHER SOURCES OF ORBITAL TORQUE 

Here we consider other perturbations that could induce 
changes in S2's orbital angular momentum, potentially com- 
plicating the signal from an IBH. We find that almost all 
such alternative perturbations are small compared to the 
torque produced by an IBH. 

Angles quoted in this section are intrinsic, not astro- 
metric. As noted above, current data are able to determine 
the orbital angles (tu, fi, i) of S2 with an accuracy of about 
one degree. Changes induced by an IBH per orbit of S2 are 
10"^deg < A(i, Q) < Ideg (Figure 2). 

Spin of the MBH. Frame-dragging effects from a spin- 
ning MBH include an additional in-plane precession term as 
well as a precession of the orbital plane. Defining i' and fi' 
as the inclination and nodal angle of S2's orbit with respect 
to the MBH's equatorial plane, to lowest FN order, frame 
dragging induces changes 



Auj = — 2yljcosi 
AO.' = Aj 



(14a) 
(14b) 



per revolution in the angle of periapse and the line of nodes, 
respectively, where 

-| 3/2 

(15a) 



A, = 



4tvx 

r-3 



GA/mbh 



(1 -e2)a 
0.115' (1 



2\-3/2 / a 

e ) X 

' \^ mpc 



-3/2 



(15b) 



and Y 1 is the dimensionless spin of the MBH (e.g. 
iMerritt et aLll2010l ). The orbital inclination i' remains un- 
changed. For S2, the spin contribution to advance of the 



periapse is ~ 1% of the Schwarzschild contribution even 
for X — 1 S'lid so can be ignored. The nodal advance is 
AQ' ~ 0.002 xdeg, too small to be detectable in the next few 
decades of monitoring, and smaller than the changes induced 
by an IBH in almost all of the runs carried out here (Fig. 2). 
Effects of frame dragging are only likel y to be important fo r 
stars much closer to SgrA* than S2 I Merritt et al.|[2010l ). 
(Frame dragging could nevertheless be relevant to the orbit 
of an IBH. For the IBH orbit with smallest a (0.3 mpc) and 
largest e (0.9) considered here, the precession time drops to 
~ SOOx"' yr.) 

A stellar bar. If the gravitational potential due to the 
distributed mass is appreciably non-spherical, the resultant 
torques could affect all of the orbital elements of S 2 aside 
from a. For instance, the nuclear bar described by lAlardI 
(|200lh has been modelled as a triaxial spheroid with central 
densi ty ~ 15OM0pc~'' ([Rodrigucz-Fcrrrandez & Combei 
[200^ . In the most extreme case, the nuclear star cluster 
(NSC) of the Milky Way could be stratified on triaxial ellip- 
soids at all radii. A homogeneous, non-rotating triaxial bar 
induces changes in the inclination and nodal angle of a test 
star orbiting near the MBH, with characteristic time scale 



(16) 



Ti - T, PGpb 

(*Mer ritt Sz Vasilievll2O10l . equations 11-15). Here, pt is the 
density of the triaxial component, and {Tx,Ty,T^) are the di- 
mensionless coefficients, of order unity, that defi ne the shape 
of the triaxial component l|Chandrasekhatll 19691 ) ; the torque 
in the principal plane is proportional to Ti — Tj , etc. For 
S2, the angular reorientation over one orbit due to torques 
from a triaxial bar would be of order 



Acp ^ 0.2mas {T, - Tj\ 



Pb 



(17) 



100 Mqpc-3^ 

the approximate 



undetectable even if pt ~ 10^ MqPc" 
density of the NSC at 1 pc from SgrA*. 

Resonant relaxation: Discreteness in the distribution 
of stars and stellar remnants is also a potential source of 
torque. "Vector resonant relaxation" produces changes in 
the direction of the orbital angular momentum of order 

^/Nm 



A0: 



(18) 



Mmbh 

per radial period, where m is the mass of a background star 
and N is the number of such stars within S2's orbit. Writing 
A'^ — M^/m, with M* the distributed mass within S2's orbit, 
this is 

1/2 



A(j) 



0.003 deg ' 



(19) 



1/2 



-1/2 



10-2 J V4xl06m/ ^'^^^ 

Since < IO-^Mmbh (iGillessen et al.ll2009l ). this is smaller 
even than the change due to frame-dragging unless the mean 
perturber mass m S> Mq . 



6 DISCUSSION AND CONCLUSIONS 

We have shown that an IBH orbiting Sgr A* at the distance 
of the S-stars can cause observable deviations in the orbit of 
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Figure 11. Root mean square velocity of the MBH as a function 
of the binary mass ratio. The dotted line repr esents the 3o" lim it 
on the 3D proper motion of SgrA* derived from lReid et al.l ||2009|) . 



star S2. The perturbations potentially affect all the orbital 
elements, but the key signature would be a change in S2's 
angular momentum (eccentricity, orbital plane) due to the 
non-spherically-symmetric forces from the IBH. In particu- 
lar, we find that an IBH more massive than ~ 1000 Mq at 
a distance of ~ 1 — 5 mpc is potentially discoverable at the 
next periapse passage of S2. 

The presence of an IBH companion to the MBH might 
produce an observable feature in the radio observations of 
SgrA*. In particular, some combinations of the binary pa- 
rameters considered here might produce a reflex motion of 
the MBH which is larger than the currently available limits 
on the proper motion of SgrA*. Such combinat ions of pa- 
ramet ers could therefore be considered unlikely. iReid et aU 
lj2009l ) report a proper motion of (7.2 ± 8.5) kms"'^ in the 
plane of the Galaxy and of (— 0.4 ± 0.9) kms~^ in the di- 
rection perpendicular to the plane. Since it appears unlikely 
that the motion of the MBH lies primarily in the Galactic 
plane, we adopt the value of the perpendicular component 
of the velocity as our fiducial value. In Figure[Tl]we compare 
the 3D root mean square velocity of the MBH obtained from 
the simulations with the observational 3a upper limit. We 
find that the motion of the MBH induced by the orbiting 
IBH is generally smaller than the 3a limit derived from radio 
observations of SgrA*. Only for q > 10"'' and a — 0.3 mpc 
does the simulated motion exceed the limit. Of course, if 
the motion induced by the IBH were to lie primarily in the 
Galactic plane, it would be very hard to detect via radio 
observations, even for large values of q. 
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